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Abstract 

We find evidence that a certain class of reaction-diffusion systems can exhibit chemical turbulence 
equivalent to Nikolaevskii turbulence. The distinctive characteristic of this type of turbulence is 
that it results from the interaction of weakly stable long-wavelength modes and unstable short- 
wavelength modes. We indirectly study this class of reaction-diffusion systems by considering an 
extended complex Ginzburg-Landau (CGL) equation that was previously derived from this class 
of reaction-diffusion systems. First, we show numerically that the power spectrum of this CGL 
equation in a particular regime is qualitatively quite similar to that of the Nikolaevskii equation. 
Then, we demonstrate that the Nikolaevskii equation can in fact be obtained from this CGL 
equation through a phase reduction procedure applied in the neighborhood of a co dimension- two 
Turing-Benjamin- Feir point. 
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The onset of spatiotemporal chaos is an important subject in the study of dissipative 
systems Several yeas ago, a new mechanism causing the onset of spatiotemporal 

chaos was discovered by Tribelsky et al. \^ for the Nikolaevskii equation, 



d t u = -d 2 [e - (1 + d 2 r ) 2 ]u - {d r uf. (1) 
This equation was originally proposed to describe the propagation of longitudinal seismic 



waves in viscoelastie media p. Its uniform steady state // is unstable with respect 
to finite-wavelength perturbations when the small parameter e is positive. However, this 
instability does not lead to spatially periodic steady states, because the equation possesses a 
Goldstone mode, due to its invariance under transformations of the form u — > u + const., and 
the corresponding weakly stable long-wavelength modes interact with the unstable short- 
wavelength modes. As a consequence, spatially periodic steady states do not occur, and 
instead spatiotemporal chaos occurs supercritically. This chaos is called 'Nikolaevskii tur- 
bulence'. Its properties have been investigated by several authors 0. Ifl lltl llll Il2 . . 

Although it is conjectured that spatiotemporal chaos exhibiting a similar onset appears in 
various systems, experimentally only one such phenomenon has been observed to this time, 



complex electrohydrodynamic convection (also called 'soft-mode turbu 



ence'), discovered in 



homeotropically aligned nematic liquidcrystals by Kai et al. Q 0, 0, HQ. 

Similar 

onset has also been studied numerically in systems exhibiting Rayleigh-Benard convection 
under free-free boundary conditions by Xi et al. ^| , and the possibility for the existence of 
this type of turbulence in reaction-diffusion systems has been investigated by Fujisaka and 
Yamada (^(J and independently by Tanaka and Kuramoto j^]. 

In this letter, we present further evidence for the ubiquity of the type of spatiotemporal 
chaos described above. We find evidence that the chaos exhibited in a particular regime by 
a complex Ginzburg-Landau equation with nonlocal coupling j^jj, called a nonlocal CGL 
equation, is equivalent to Nikolaevskii turbulence. This suggests that a certain class of os- 



cillatory react ion- diffusion systems can also exhibit this type of chaos 



Decause the nonlocal 
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221. These reaction- 



CGL is a reduced form of this class of reaction-diffusion systems 
diffusion systems are such that the chemical component constituing local oscillators are only 
weakly diffusive, while there is an extra diffusive component introducing an effectively non- 
local coupling between the oscillators. Such a situation, in which the coupling between the 
local oscillators is mediated by a diffusive substance, has been observed in some systems 
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studied experimentally: biological populations, such as cellular slime molds [23j and oscil 



n 

lating yeast cells under glycolysis 24], catalytic CO oxidation on metal surfaces 



3, and 

the Belousov-Zhabotinsky reaction exhibited by a system dispersed in a water-in-oil aerosol 
OT microemulsion [26]. Because these systems possess such coupling, it is conjectured that 
under certain circumstances they could exhibit spatiotemporal chaos in the same class as 
Nikolaevskii turbulence. 

Our starting point is the following nonlocal CGL equation with complex amplitude A: 

d t A = A - (1 + ic 2 )\A\ 2 A + (5 1 + i5 2 )V 2 A 

+K{l+i Cl ) Jdr'G(r - r')[A(r') - A(r)}. (2) 

Here G is a coupling function, and c\, c 2 , Si, 5 2 and K are real parameters. In a previous pa- 



per 



2lj , this equation was derived as a generic reduced form of the class of reaction-diffusion 
systems discussed above. A simple example that belongs to this class is the following set of 
equations, describing a hypothetical extended Brusselator: 

d t X = a-(b+ l)X + X 2 Y + DV 2 X + kc x S, (3) 
d t Y = bX - X 2 Y + DV 2 Y + kc Y S, (4) 
rd t S = -a-S + V 2 S + X, (5) 

where the field variables X and Y represent limit-cycle oscillators existing immediately 
above the Hopf bifurcation, and the additional chemical component S introduces an effective 
nonlocal coupling between these oscillators. The strength and anisotropy of this coupling 
are represented by the parameters k and cx,y respectively. The Hopf bifurcation parameter 
fi is written in terms of the parameters a and b as = (6 — b c )/b c , where b c = 1 + a 2 . The 
coefficient D is the diffusion constant for X and Y, and the quantity r is the time constant 
of S. When these couplings are as weak as the oscillation, i.e., when fi ~ 0(D) ~ 0(k), this 
system can be reduced to Eq. (J2J). 

Equation (j2J) is invariant under transformations of the form A — > Ae lc , with real constant 
c, which implies that a uniform mode is neutrally stable. This is a result of the spontaneous 
breaking of time translational symmetry corresponding to the Hopf bifurcation in the original 
react ion- diffusion systems. Also, the uniform oscillating solution of Eq. (J2J possesses a 
Turing instability |2j| in a certain parameter region, as shown in Fig. ^ in contrast to the 
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FIG. 1: Parameter space (K, Kc\) in the neighborhood of a co dimension- two Turing-Benjamin- 
Feir point. The uniform ocsillating solution of the nonlocal CGL is stable in domain S, possesses 
the Turing instability above the curve It and the Benjamin- Feir instablility above the line Ibf- 
The parameter values are c2 = 1, 5i = 0.3, 62 = 0, and G(r) = p/2 e~ /9 ' r ', with p = l + 0.905i. The 
inset plots the stability eigenvalue as a function of the perturbation wavenumber for the uniform 
ocsillating solution at the point (K,Kc\) = (0.05,-0.4345), indicated by the cross. The largest 
eigenvalue is nearly equal to 0.0005 at the Turing wavenumber, which is approximately 0.6. This 
dispersion curve is qualitatively the same as that of the Nikolaevskii equation. 



ordinary diffusive complex Ginzburg-Landau equation, which possesses only a Benjamin- 
Feir instability. Thus, Eq. (j2j) in this parameter region is characterized by weakly stable 
long- wavelength modes and unstable short- wavelength modes, like Eq. (JI|) . Therefore, it 
has been h ypo thesized that Eq. (j2J exhibits spatiotemporal chaos similar to Nikolaevskii 



turbulence 21]. We confirm this hypothesis in the following. 

For simplicity, we fix the parameter values as C2 = 1 and 62 = and consider only 
the case of one spatial dimension. In this case, the coupling function G{r) is given by 
G(r) = p/2 e _p ' r ', with p = 1 + ir) and 1 > 77 > 0, where the form of G(r) is derived from the 
original reaction-diffusion systems [21]. However, the following analysis is applicable also in 
higher-dimensional situations, in which G takes other forms. 

A typical spatiotemporal pattern exhibited by Eq. (J2J in the slightly Turing-unstable 
regime close to the Benjamin- Feir criticality is shown in Fig.0 where we see long-wavelength 
modulation of the Turing pattern. The corresponding spatial power spectrum is displayed in 
Fig. 01 This spectrum is found to have characteristic peaks at the Turing wavenumber and 
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FIG. 2: Gray level plots of \P(r, t) = d r axg[A(r, £)], which corresponds to d r u of the Nikolaevskii 
equation Eq. Q). The width of the longitudinal stripes corresponds to the Turing wavelength. The 
parameter values are indicated by the cross in Fig. ^ The configuration displayed here represents 
the result obtained after evolution from an initial configuration (which is the uniform ocsillating 
solution with a slight perturbation) over a time interval of length 10 5 in a periodic one-dimensional 
system consisting of 2 10 spatial points separated by a distance of 0.3. For the numerical scheme, 
we used a pseudo-spectral method with an explicit fourth-order Runge-Kutta scheme of time step 



its harmonics. This feature is also seen in the spectrum of Nikolaevskii turbulence, shown 
in the inset of Fig. El This suggests that in the regime we consider, the spatiotemporal 
chaos exhibited by Eq. (j2J) is in the same class as Nikolaevskii turbulence. In the following, 
we show that, in fact, Eq. can be obtained from Eq. (J2J) by means of a phase reduction 
technique 

In the phase reduction, the diffusion term and the nonlocal-coupling term of Eq. (J2J) are 
treated as perturbations. Using Floquet theory, we can calculate the operator Z(0) that 
projects the perturbation onto the limit cycle A ((f)) = e - ^ of the unperturbed system: 
Z((f>) = (cos — sin 0, — sin — cos 0). Hence, the phase dynamics of Eq. (J2J) obey 



0.01. 



d t <p 



l + 5i<^0-5i(<9 r 0) 2 




(6) 
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e ' r '[sin(0 — r] \r\ + arga) 
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sin(— rj \r \ + argo;)] 



(7) 
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FIG. 3: Spatial power spectrum of ^f(r, t). This spectrum was calculated using a system size 
(0.3 x 2 12 spatial points) larger than that of Fig. [21 and averaged over a long time, excluding 
an initial transient of length 5 x 10 4 . The inset displays typical spatial power spectra obtained 
from the Nikolaevskii equation, Eq. with e = 0.02 (bottom) and the well-known Kuramoto- 
Sivashinsky equation (top), for the sake of comparison. The peaks are characteristic of Nikolaevskii 
turbulence. The spectrum of the nonlocal CGL exhibits similar peaks at the Turing wavenumber 
and its harmonics. 

where 0(r, t) and 0(r', t) are abbreviated as and 0', respectively, and a = [—c\ — 1 +z(l — 
Ci)]p. In the integrand r(0— 0', r—r') there is a dependence on the quantity 0— (0'+7/|r— r'|); 
i.e., the time evolution of the phase at r depends not on the difference between the phases at 
r and r' but on the difference between the phase at r and the phase at r' plus the quantity 
r]\r — r'\. This might seem paradoxical for the following reason. If the spatial interactions 
are weak, then d t (j) — 1, and hence 0(r',t) +r]\r — r'\ ~ 0(r', t + 1 ) with t = r]\r — r'\. Thus 
the phase at r interacts with the future phase <p(r',t + t ), where to ^ because rj > 0. To 
solve this paradox, note that the original react ion- diffusion fields oscillate roughly as cos $ 
with $ = uot — /i0, where the Hopf frequency uq is a finite positive value, and the Hopf 



bifurcation parameter /i is an extremely small positive value [2l|. Thus, $(r, t) interacts 
with $(r', t — ti) with ti = /ztoA^o if we ignore terms of 0(fi 2 ). Hence the field at r interacts 
with the past field at r', as ti > 0. In this sense, rj producing the imaginary part in the 
coupling function G in Eq. (J2J) causes the phase coupling to be implicitly delayed, with the 
delay proportional to the distance between interacting oscillators. 

Now we rescale the phase equation derived above. First, we introduce a variable ip defined 
as ip = — t. The evolution equation for ip is the same as that of 0, except for the absence 
of the term 1 on the right-hand side. Instead of Eq. we refer to the equation for if) 
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as the phase equation in the following. In the long- wavelength limit, where the validity of 
the phase description is ensured because the amplitude-like modes can be safely ignored, 
the phase equation can be expanded in a power series in d 2 , owing to the symmetry of the 
system with respect to reflection through r = 0: d t ip = (A 2 c^ + A 4 <9^ + X 6 df + ■ ■ ■ )tb + N. 
iere N represents the nonlinear terms, and A n are constants derived from c\, 8\, K and r) 
. We consider the Turing-unstable regime close to the Benjamin- Feir criticality, where 
A 2 , A 4 , A 6 > 0, A 2 ~ 0(v 2 ), A 4 ~ and A 6 ~ 0(1), with the scaling parameter v = +0, 

and 4A 2 A 6 = (1 — e)A 4 with a positive constant e C 1. Under these conditions, the higher- 
derivative linear terms, d™ib with n > 8, are much smaller than the other linear terms and 
can be ignored because the characteristic spatial scale of ib is ~ 0{y~ x l 2 \ Furthermore, 
because ib itself has a characteristic small magnitude depending on v, it is reasonable to 
assume that the largest nonlinear term in N is '-/(drip) 2 , where 7 is a constant derived from 
Ci, Si, K and i] (^j]. In fact, when we rescale r, t and ip as r ~~ > ? — yj ^ r -> t ^ t = and 
ib — > tA = ~ 4 % Xe , ib, we can obtain the following scale- free equation from the phase equation: 

d^ = -dl[e-{l + dlf]^-{d^f. (8) 

This is identical to Eq. (0). Here, note that ib satisfies the scaling relation 

ip(r,t) = v 2 i){y xl2 r,vh). (9) 

This spatiotemporal scaling of the phase is completely different from the scaling relation in 
the case that the Kuramoto-Sivashinsky equation is derived, ib(r,t) = ^{^ l ^ 2 T,^ 2 t), where 
£ is a parameter that represents the Benjamin- Feir criticality 

The above parameter conditions resulting in the derivation of Eq. (JHj) can be written in 
terms of the parameters of Eq. (J2J) as follows. Using ( = arctanr/, where 7r/4 > £ > 
because 1 > r\ > 0, Eq. (JHJ) is derived from Eq. (J2J under the following conditions, with a 
sufficiently small positive constant v when ( 7^ 0,7r/8: 

K(l + i Cl ) = Texp[i(4C + 37r/4-i/)], (10) 

where 

Si 



T 



v^cosC) 2 



. . , (1 — e)(sin vY 
sin(2C - v) + V 



4 sin(2C + u) 

Here, the parameter e is the same as that in Eq. (jHJ), and the critical value of T is T 
5 1 [v / 2(cosC) 2 sin(2C)]" 1 . 



(11) 



In order to observe the same spatiotemporal chaos as that of Nikolaevskii turbulence in 
our reaction-diffusion systems, we can analytically tune these parameters to satisfy the above 
conditions. In particular, the time constant r of the additional chemical S (e.g., see Eq. Q) 
should be nonzero, because ( = [arctan(o; r)]/2 7^ 0, where Uq is the Hopf frequency [21]. 
This reflects the importance of the imaginary part in the coupling function G of Eq. (J2J), i.e., 
the implicit delay of the phase coupling discussed above. This contrasts with the situation 
found in previous studies of similar reaction-diffusion systems, in which the limit r — ► was 



taken and S was eliminated adiabatically 
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28, 




29, 30L 13 li . In addition to the condition 



r ^ 0, we should also have r 7^ 1/&>q, because ( 7^ ir/8. This implies that the charcteristic 
time of the field S is not comparable to that of the local oscillator, (X,Y). 

In conclusion, we have found evidence that the nonlocal CGL and the corresponding 
class of reaction-diffusion systems in a certain regime exhibit turbulence that is equivalent 
to Nikolaevskii turbulence. This result supports the conjecture of the ubiquitous nature 
of spatiotemporal chaos caused by the interaction between weakly stable long-wavelength 
modes and unstable short-wavelength modes. We have confirmed through numerical cal- 
culation that these chaotic states are structurally stable for the nonlocal CGL not just at 
the codimension-two Turing-Benjamin- Feir point but also in a finite neighborhood around 
it, at points where the CGL equation does not reduce exactly to the Nikolaevskii equa- 
tion. Furthermore, we believe that similar spatiotemporal chaos would be found even if 
the phase-shift invariance were sligthly broken in the reaction-diffusion systems; i.e, this 
chaos should persist even if the Goldstone mode were lost. This conjecture is based on the 
observation that soft-mode turbulence exists even if one applies a small magnetic field that 
slightly breaks the arbitrariness of the azimuthal angle of directors bent through a Freed- 
ericksz transition in homeotropically aligned liquidcrystals 3^. Finally, we note that the 
phase equation derived from the nonlocal CGL is a useful model for studying quantitative 
features of Nikolaevskii turbulence that has not yet been investigated sufficiently Because 
this phase equation covers the transition region between Nikolaevskii turbulence and the 
well-known Kuramoto-Sivashinsky turbulence, it can be used not only for comparing these 
two types of turbulence but also for exploring the transition between them. This should lead 
to deeper understanding of Nikolaevskii turbulence and also a broader class of turbulence 
caused by interactions among modes with vastly different length scales. 
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